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Recent data allow a new parametrization of the neutron charge form factor GJ. A parameter- free 
quark-model relation between and the N— >A quadrupole form factor Gq^* a is used to predict 
Gc2~* A from data. In particular, <r^> is related to iV^A quadrupole moment Qjv^a, while 
<r*> connects to the N—*A quadrupole transition radius <r%^ A >. From the latter we derive an 
experimental value for the charge radius of the light constituent quarks r lq — 0.8 fm. Finally, the 
C2/M1 ratio in pion electroproduction is predicted from the elastic neutron form factor data. 
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Intrinsic nucleon structure - The nucleon is a complicated many-particle system composed of valence quarks, which 
carry the quantum numbers, and nonvalence quark degrees of freedom, which describe the cloud of quark-antiquark 
(qq) pairs and gluons. The constituent quark model (CQM) with two-body exchange currents describes both these 
aspects of nucleon structure jlj . One-body currents describe the interaction of the photon with one valence quark at 
a time. Two-body exchange currents are connected with the exchange particles (gluons, pions) and with qq pairs (see 
Fig. |l|(b-d)). Nucleon properties which are dominated by two-body exchange currents show their common dynamical 
origin in analytical interrelations. 

In this paper, a quark model relation between the neutron charge form factor and the quadrupole transition 
form factor Gg 2 ~^ A is used together with a parametrization of new G n E data to predict G^ A and the C2/M1 ratio 
usually measured by pion electroproduction. The spin-isospin dependence of the two-body charge operator p[ 2 ] , e.g. 
for the gluon, can be written schematically as 



(Ti-CT^iq) - ^ jjtr i x <x,] 2 x F 2 (q)] 



(1) 
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where q is the three-momentum transfer, ej and <Xj the quark charge and spin. Explicit expressions can be found in 
Ref. fl. 
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FIG. 1. Feynman diagrams of the four- vector current J M with photon coupling to (a) one-body current J?,, and to (b-d) 
two-body gluon and pion exchange currents . 



As a result of (i) the dominance of pp] IIJ4] and of (ii) the spin-isospin structure in pp] and in the N and A wave 
functions, a connection between the neutron charge form factor G^(q 2 ) and the N^A quadrupole transition form 
factor Gc 2 ~^ A (q 2 ) emerges (U as follows: 

G£r A (q 2 H-^G'j;( q 2 ). (2) 

The derivation is independent of the spatial part of the quark wave functions and holds irrespective of whether gluon 
or pion exchange or both are employed. One-body contributions to G^ and Gq 2 ~* a arise only through excited-state 
admixtures which amount to less than 20% of the empirical values if a quark core radius consistent with the excitation 
spectrum is used. Three-body corrections are estimated to contribute less than 30% of the two-body currents 
using a QCD parametrization Q and a large 1/N C approach ||. These approaches also show the dominance of the 
two-body exchange currents for both observables. 

The dominance of nonvalence quark degrees of freedom in GJ and Gq 2 ^ A is not specific to this quark model. Also 
chiral approaches, e.g. the Skyrme model R], the er-model H, the Nambu-Jona-Lasinio model M, and the chiral 
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soliton model with quarks jlO| predict that the neutron charge radius <r%> and Qjv->Aj the transition quadrupole 
moment, are dominated by nonvalence quark degrees of freedom. In these models, the valence quark contribution 
to <r 2 > and Qn^a is comparatively small. The quark model, the chiral approaches fLl| , and chiral perturbation 
theory [fl2f have the same underlying group theoretical structure, i.e. spin-isospin symmetry. Eq.(^) has been derived 
within the constituent quark model for the case of dominance of two-body operators. However, a larger range of 
validity is conjectured because of stringent constraints on allowed operator structures in Eq.(Q) due to the scalar 
nature of the charge operator and the spin-isospin symmetry of the N and A wave functions. 

For further interpretation of the new parametrization we rely on explicit expressions for the lowest moments of 
GJ. and G^^ A calculated in the quark model with two-body currents. In this model the neutron charge radius 
can be expressed in terms of the empirical N-A mass splitting and the quark core radius b 



<rl> = -b' 



Ma — M, 



N 



M. 



N 



(3) 



Inserting the experimental neutron charge radius <r^>=-0.113 fm 2 and the experimental N-A mass splitting 
into Eq.(|) one obtains b = 0.60 fm for the quark core radius, which measures the spatial extent of the valence quark 
wave function. A related quantity of interest is the slope of Gq^ a at q 2 =0, which is given by [^|,[| 



N—>A 



>= — b 
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where r 2 is the charge radius of the constituent quark. In the same model we find G 



>A (q 2 ) 



-V2 G^(q 2 ) 



Moments of form factors 
into a Taylor series 



First, at low momentum transfers Q =—q any form factor, (here G^), can be expanded 



G|(Q 2 ) = G1(0) 



dG£ 



dQ 2 
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Second, in the Breit-frame (Q 2 =q 2 ) the form factor G^. is related to the spatial charge distribution p(r) through the 
Fourier transform as 



G'KQ 2 ^ 



p(r)e^ qr dr = 4tt / pM ^^r 2 dr 

Jo q r 



(6) 



A series expansion of sin(qr) leads to the well known relation between the derivatives of the electric form factor and 
the radial moments of the charge distribution 
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The r.h.s. are obtained by taking the first and second derivative of the common parametrization of G^ 



G|(Q 2 ) 



1 + dr 



G_d(Q 2 ) 



(9) 



(10) 



with G D (Q 2 ) = 1/(1+Q 2 /A 2 ) 2 , A 2 =0.71(GeV/c) 2 , and r=Q 2 /(4M 2 ) (see refs. @,[|), and with p n the neutron 
magnetic moment. This simple parametrization guarantees the proper behaviour of G^ at Q 2 =0 (zero net charge) 
and at Q 2 — >oo (quark counting rules). 

Similarly, the low momentum expansion of Gq 2 ~^ A can be written as 



c£r A (Q 2 ) 



lN->A 



< r 



N — >A 



> 



Q 2 



(ii) 



2 



Referring to the connection between Gq 2 ~* a and (Eq.|J) and equating the respective coefficients of Eqs.([Tl|) and (||) 
one finds with the help of Eqs.(0^) 

<r 2 n > 3a /x„ 

Qn - a =^r = (12) 

<r 2 > _ 6 <vj> _ 3 <r*> _ 




Q 2 (GeV/c) 2 



FIG. 2. Recent measurements of G^ are given by the full points (see Table H). The Platchkov parametrization [Jl5| (dashed) 
is shown together with a fit to the new data without (dot-dashed) and with (full line) inclusion of <r^> from Ref. Ea|. The 
insert shows the neutron charge distributions resulting from the new (full line) and the Platchkov (dashed line) parametrizations. 



New GTg data and fits - There has recently been increased interest in the neutron electric form factor itself as 
a testing ground for nucleon models. Furthermore, G^ must be known accurately in order to analyze for example 



parity violation and N— +A transition experiments. The previous evaluation by Platchkov et al. |15| of elastic electron- 
deuteron scattering data is model-dependent because the analysis relies heavily on the deuteron wave function when 
subtracting all magnetic form factors and the charge contribution of the proton in order to isolate G^ from the 
measured cross sections. The analysis employing the Paris potential was then considered as the most reliable p5| , 
other choices for the NN potential yielded values for G^ differing up to 100%. The results were parametrised using 
Eq.(|l0[) with the two parameters fitted to a=1.25 and d=18.3. An earlier analysis obtained the values a=1.0 and 
d=10.7, respectively. 

TABLE I. GJ data resulting from recent experiments are selected and are listed with statistical (A s t a t) and total 
(Atot=A s tat + A sys t) errors. The <r„> values from electron-neutron scattering are given in the last line. 



Q 2 a ) 




A s tat 


Atot 


ref. 


0.15 
0.21 
0.34 
0.67 


0.0480 
0.0660 
0.0611 
0.0520 


0.0065 
0.0150 
0.0069 
0.0110 


0.0118 
0.0190 
0.0129 
0.0160 
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Astat b ) 


At„t b ) 






0.0 


-0.113 


0.0026 


0.0060 
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a ) in units of (GeV/c) 2 b ) in units of fm 2 
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In modern double-polarisation experiments, where the ratio of G^/G^ is measured in quasifree kinematics, the 
NN-potential dependence is greatly reduced. The remaining uncertainty resides mainly in the corrections for final 
state interaction (FSI) and to a lesser extent for meson exchange currents (MEC) between the nucleons. Recent 
experiments have been performed at Bates, Amsterdam and Mainz by using polarised electrons scattering either 
from unpolarised deuterium and measuring the recoil polarisation of the neutron |l6|-|il^] , or by measuring helicity 
dependences on polarised 3 He po|-^3|. Only four of these data with sufficient statistical and systematic accuracy 
have been selected (Table ||) and are displayed in Fig. ^. A further selection criterium was the coincident detection of 
neutrons and scattered electrons which enables a check on the quasifree scattering mechanism. 

The three deuteron data [p~6| — p~S|| rely on corrections based on calculations by Arenhovel [ p4| . Corrections for the 
3 He data at Q 2 =0.35 (GeV/cp ^,^| are expected to be quite large. Therefore, these 3 He data were not considered 
in the present fits. Note however, that the correction at Q 2 =0.67 (GeV/c) 2 is estimated to amount to only 10 %, 
which has not been applied to the data but is included in the systematical error p(i| ]. 

Statistical errors are given in column 3 of Table [I]. The sum A to t=A sta t+A syst of statistical and systematical errors 
is shown separately in column 4. In view of the discussion above, we will use A to t in this analysis. The data are shown 
in Fig. H together with the (dashed) curve based on the Platchkov parametrization with the Paris potential (ill] . We 
note that the new data lie significantly above the previous evaluation. 

Parameters a and d of Eq.(|l^) have been obtained by a "downhill simplex" fit using the four data points from 
quasifree scattering and the < r 2 >-value ]l3| obtained from thermal neutron electron scattering. The results and 
their standard deviations due to the fit with A tot are compiled in Table O, together with the total x 2 - For reference, 
the parametrizations of Galster [ ft4[_a nd Platchkov are listed in the last two lines of Table || showing large values 
for x 2 ■ The first two lines of Table III present the fitting results for different conditions as marked in the first column. 
From Eq.fj?]) it is evident that parameter a is determined exclusively by <r 2 >, which leads to a strong reduction of 
Act when <r 2 > is included in the fit. The result including this constraint (line 2 of Table ||) is regarded as the most 
reliable, considering the present status of the data base. It is represented by the full curve in Fig. ^. The influence 
of omitting the < r 2 >-datum in the fit (line 1 of Table |l|) is shown by the dot-dashed line. From this fit a neutron 
charge radius of < r 2 >=-0. 178(27) fm 2 would be deduced which is ruled out by the analysis of Ref. JT^I - Despite its 
large uncertainty the parameter d is much smaller than previously assumed ]15[ which is reflected by smaller values 
for G^/'(0) and 1Z as shown in the last two columns. 

TABLE II. Parameters a and d from fits to the data of Table [j] together with \ 2 are shown in the first two lines; values for 
Gb"(0) (Eq. ^) and the ratio 1Z (Eq. ^) are also given. The next three lines contain results where one additional ficticious data 
point with At ot =10 % at Q 2 =0.9 (GeV/c) 2 is included in the data set for which the assumed G^-value is given in col. f . For 
reference the results using the parameters of Refs. [B and |f5| are shown in the last two lines. 



comments 


a ± Aa 


d ± Ad 


x 2 


G n E " a ) 


TZ b ) 


c ) 






1.415 (0.213) 


9.06 (3.30) 


0.09 


-8.26 


-10.77 








0.898 (0.044) 


2.74 (1.99) 


0.63 


-3.50 


-7.19 


G n E = 


0.04 


0.903 (0.040) 


4.30 (0.77) 


1.30 


-3.95 


-8.07 


G n E = 


0.05 


0.898 (0.041) 


2.78 (0.66) 


0.63 


-3.51 


-7.21 


/"in 


0.06 


0.894 (0.040) 


1.75 (0.61) 


1.07 


-3.21 


-6.63 


ref. [ 
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1.00 ( - ) 


10.70 ( - ) 


10.53 


-6.35 


-11.69 


ref. [ 


L5 




1.25 ( - ) 


18.30 ( - ) 


62.65 


-10.83 


-16.00 



a ) in units of (GeV/c) 4 c ) <r^> datum omitted 

b ) in units of (GeV/c)- 2 
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With the present fit the inverse Fourier transformation of E g . <S ) leads to a neutron charge distribution (full line in 
insert of Fig. ||) quite different from the previous results (Ref. flEfl, dashed line), and thus to different moments. The 
inner region of positive contributions is compressed and the zero crossing point is shifted by about 0.2 fm towards the 
center, compared to the old distribution. Note, that the new zero coincides with the value for the quark core radius 
b derived from Eq. (||). 

A better determination of parameter d, which is crucial for the second derivative of , calls for additional data 
at high Q 2 where corrections for FSI and MEC are less important. Although the cross sections decrease with Q 2 , 
a statistical error A stat =5 % is attainable and the chosen error A tot =10 % seems to be realistic. Simple values for 
G^ have been assumed to demonstrate possible variations of d. In Table || we present three more fits each with one 
fictitious data point added at Q 2 =0.9 (GeV/c) 2 . The additional datum at Q 2 =0.9 (GeV/c) 2 reduces Ad by a factor 
3. Because the highest data point at Q 2 =0.6 (GeV/c) 2 does not include any FSI correction, a further decrease of 
parameter d towards a value of 2 seems likely. 

C2/M1 : — ' — ' — ' — i — 1 — 1 — ' — i — ' — ' — 1 — i — ' — 1 — 1 — i — 1 — T 

o.oo ; ; 




FIG. 3. Ratio C2/MX from the present fit to (solid curve) and from calculations of Ref. |3J (dashed curve) in comparison 
with experimental results taken from Refs. 



Results - The fits of Gj| and in particular the first extraction of the fourth moment <r„>=-0.32(8) fm 4 permit a 
prediction for <r 2 v ^ A >=0. 84(21) fm 2 where the error is determined only through Ad (see Eq.@); with the additional 
sixth data point the error of <r%^ A > would be reduced to 0.07 fm 2 . The variation of d and thus 1Z with the sixth 
data point is indicative of a possible 10 % change of <rj^^ A >. By use of Eq.(H) and a quark core radius 6=0.6 fm 
we obtain for the charge radius of the light constituent quarks of r 2 ? =0.64 fm 2 . This value is somewhat larger than 
the value derived from the vector dominance model, which would give r 2 g =6/m 2 =0.4 fm 2 . 

Finally, in the quark model with two-body exchange currents the C2/M1 ratio can be expressed as a ratio of the 
elastic neutron charge and magnetic form factors || 

<X =M GgT A (q 2 ) _ M N G|(q 2 ) 

Ml N 6 G$r A (q 2 ) 2^ G^(q 2 ) [ > 

Using G^ f =/i rl Gz) we compare the predictions of Eq.(|l4|) with the direct measurements p5| , p6| of the C2/M1 ratio 
in pion electroproduction. Starting from the new parametrization of G^ we find the ratio in the range of -0.03 to 



-0.08 (full line in Fig. |3|). The data |25|26|] in general are in the range of -0.046 to -0.06. The agreement is surprising, 
which hints at a general validity of the relation between G^/G^f and C2/M1 (Eq.|l4|). A constituent quark model 
calculation of this ratio (dashed line) comes lower than the parametrization. Differences between this calculation 
and the C2/M1 data may be explained by possible background amplitudes contributing to the iV^A transition. 
Other constituent quark model calculations based on one-body currents alone, e.g. Ref. [^7j, predict much smaller 
ratios, mostly in the range of -0.005 to -0.02 (see also refs. in p5|). 

In summary, we have performed a fit to the most reliable data for G^ below Q 2 <1 (GeV/c) 2 . We use the ob- 
tained parametrization in combination with relations originally derived in the constituent quark model to predict 
the quadrupole transition form factor Gq 2 "^ a and its leading moments Qn^a and <r%^ A >. Eq. @ provides a 
determination of C2/M1 through the elastic neutron form factors which agrees well with C2/M1 data from pion 
electroproduction experiments. 
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